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Abstract 

We address the problem of defining the four-potential, A®(x), on the lattice, in 
terms of the natural link variables, U fl (x). 

Different regularized definitions are shown, through non perturbative numerical 
computation, to converge towards the same continuum renormalized limit. 
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1 Introduction 



Lattice QCD does not require, in itself, any gauge-fixing in order to compute physi- 
cal quantities. However, as became increasingly clear in recent years, lattice gauge- 
fixing provides a necessary instrument in the study of quantities, like the quark and 
gluon propagators, whose behavior could be relevant to the study of the confinement 
mechanism Q - 1| . 

Gauge-fixing is also a necessary ingredient in some non-perturbative renormaliza- 
tion schemes H |7| and it has been shown to facilitate the construction of composite 
fermion operators with correct chiral behavior ||, essential for the study of hadronic 
weak interaction phenomenology || - p[ . 

Once gauge-fixing has been performed, Green's functions containing the funda- 
mental quark and gluon field insertions, become accessible to non-perturbative study. 
However, contrary to what happens for quarks, a natural definition for the gluon field is 
missing. In the general field theoretical framework, as checked in perturbation theory, 
this is known not to be a real problem: any pair of operators differing from each other 
by irrelevant terms, i.e. formally equal up to terms of order a, will tend, to the same 
continuum operator, up to a constant. It is the purpose of this paper to show that this 
feature is also satisfied at the non-perturbative level, within the framework of lattice 
QCD. We will, in fact, show that different definitions of the gluon field, at the lattice 
level, give rise to Green's functions proportional to each other, thus guaranteeing the 
uniqueness of the renormalized continuum gluon field. 

The plan of the paper is as follows. In section B we recall some basic facts about 
lattice gauge-fixing. In section |3| we discuss the ambiguities intrinsic to the lattice 
definition of gauge potential and in section |4| we summarize our results. 



2 Lattice Gauge-Fixing 

The Landau or Coulomb lattice gauge-fixing procedures are well known and several 
numerical algorithms are available to this aim. In the standard approach [jl], ^| the 
functional: 

F[U n ) ee -T^Re Tr£J») (1) 

H x 

is minimized with respect to Vi{x). In eq.(|l]) V is the lattice spatial volume, T its time 
extension and U~-(x) = Q(x)U /1 (x)fl(x + /u)' is the compact SU(3) gauge field, gauge 
transformed by the local gauge transformation O(x). The extrema of F with respect to 
$7, correspond to configurations satisfying the gauge condition d^A^ = in discretized 
form, with the lattice An definition: 

_ (Ufijx) ~ Ul{x))traceless 

A^x) = 2 . agQ , /x-1,.,.4, (2) 

where a is the lattice spacing and go is the lattice bare coupling constant. The nu- 
merical behaviour of the gauge-fixing algorithm is usually monitored, as a function 
of the number of sweeps, by two quantities. The first is F[i7 n ] itself, and decreases 
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monotonically. The second, denoted by 9, is denned as: 



' = -1- Y^Tr [A(x)At(x)] , (3) 



where: 



A(s) = £(A»-Ag(x-A)). (4) 

The functional F[U ] is a lattice discretization of J d A x Tr (A"^), while 9 corresponds 
to the continuum quantity J d 4 x Tr (d^A^) 2 . By its very definition, 9, as a functional 
of f2, decreases (not strictly monotonically) during the gauge-fixing process, becoming 
zero when .F[i7 n ] gets constant: its value controls the fulfillment of the gauge condition. 
The main unsolved problems concerning gauge-fixing stem from the existence of 



both continuum and lattice Gribov copies[12]-[15] and the numerical noise that they 
can generate. 

Being a result of discretization effects, the lattice Gribov copies are of a quite 
different nature from those of the continuum, related to topological obstructions of the 
gauge-fixing condition d^A^ = O]. 

This confusing situation of the numerical gauge-fixing has been already stressed in 



the literature 16]. The real concern is, of course, the influence that these phenom- 
ena may have on the value of continuum observables, when computed through the 
intermediary of non gauge invariant quantities in the schemes referred to above. 



3 The Lattice Gauge Potential 

In this section we will discuss the problems related to the ambiguities in the lattice 
definition of the gauge potential. 

A natural definition of the 4-potential in terms of the links, U^, which represent the 
fundamental dynamical gluon variables, is given in eq.@. This definition is naively 
suggested by the interpretation of Un{x) as the lattice parallel transport operator and 
by its formal expression in terms of the "continuum" gauge field variables, Au(x) as: 

Ufj,(x) = exp(ig aA^(x)). (5) 

A formal expansion in powers of a, shows that eq.(||) defines A^x), up to terms 
formally of order a 2 . 

It is clear that the definition given in eq.((2|) is far from unique: it cannot be preferred 
to any other definition with analogue properties as, for instance: 

which in fact differs from the one given in eq.(§) by terms of O(a) that formally go to 
zero as a — > 0. 

^From the algorithmical point of view, however, the various definitions are not 
interchangeable. In fact let us see what happens if we fix the gauge of a thermalized 
configuration stopping the gauge-fixing sweeps when 9 < 10~ 14 and then define 9 as 
having the same functional form of 9, as in eq. (||), with A^ replaced by A^. 
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The values of 6 and 0' during the minimization of F are reported in Fig.[l], for 
a typical thermalized configuration, as functions of the lattice sweeps of the numer- 
ical gauge-fixing algorithm. As clearly seen 6 does not follow the same decreasing 
behaviour as 6: after an initial decrease, 6 goes to a constant value, many orders 
of magnitude higher than the corresponding value of 8. This fact has been already 
remarked in ref.[]l7j where it has been attributed to the large contribution of order a 
terms. 

This marked difference between the behavior of 6 and 9 seems to cast doubts on 
the lattice gauge-fixing procedure and on the corresponding continuum limit of gauge 
dependent operators. On the contrary, we will show, at the end of this section, that 
this discrepancy has a natural field theoretical explanation. 

The relation between the two lattice definitions A^{x) and A^x) is of the form: 

A'^x) = A„{x) + a 2 W^x) (7) 

where the Wu{%) is a dimension 3 operator 1 with the same quantum numbers of An(x). 
The operator W^(x) in eq.(^) is not renormalized, so that, while formally it would seem 
safe to neglect its contribution in the continuum limit, in fact we must take it carefully 
into account. We start constructing, out of W^{x), a finite, renormalized operator, 
W*{x) as: 



1 - C(go) 



W?{x) = Z w {g Q ,a^ R )(W^x) + -^A^x)) (8) 

where fi R is the mass renormalization scale of the theory and Zyy is a logarithmically 
divergent, subtraction dependent constant, while C, as a consequence of the Callan- 
Symanzik equation, can only depend on the bare coupling <?o[18|- Eqs.(0) and (§) show 
that 

A\(x) = C(g,)A,(x) + ^W?(x) (9) 
so that, up to terms truly of order a 2 : 

A\(x) = C(g )A IM (x) (10) 
For Green's functions insertions we have, therefore, in general: 



(■ ■ ■ A',(x) . . .) 
(...A^x)...) 



C(go) (11) 



We have numerically checked some consequences of eq.([ll]) by measuring on different 
lattices, whose characteristics are reported in Table ffl, the following Green functions 
for SU(3) in the Landau gauge with periodic boundary conditions: 

(A A )(t) = -^^Tr(A (x,t)Ao(y,0)) (12) 
(AA)W = ^^Tr(4 i (x, i )4 i (y,0)) (13) 

i x,y 



1 We are assuming that only one operator W^(x) is present, in order to simplify the presentation. In general 
several dimension 3 operators exist and their mixing has to be taken into account, but the conclusions remain 
unchanged. 
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(dAdA)(t) = -L I ^^Tr(^ t (x,t)a^(y,0)} (14) 

H,v x,y 

(A dA)(t) = * ££Tr<A o (x^VWy,0)) (15) 

(AdA)(t) = -l^jTr^lx^lU^O)) (16) 



using both A and A , as denned in eqs.(H) and (g). In eqs.(0)-(16) the trace is over 
the color indexes, \i and v run from 1 to 4 and i from 1 to 3. Here and in the following 
we define A^(t) = £ x ^(x,i). 



The correlators defined in eqs. (|12[)-(|T6|) , and in particular (AiAi)(t), are relevant to 
the investigation of the QCD gluon sector In this paper we will not be concerned 

with the interpretation of their t behavior: our aim is to show their independence on 
the ambiguities related to the definition of the gluon field. 

The proportionality factor, C(go), may depend on the direction fj,, if the lattice 
breaks cubic symmetry. In our case, as shown in Table ffl, two of the lattices have 
a time extension different from the spatial one, so that we have a coefficient Co (go) 
relating A' to Aq and a different one, Ci(go), connecting A[ to A%. 

It is worth noting that (AoAo)(t), when evaluated through A^(x), should be con- 
stant in t configuration by configuration, in virtue of the Landau gauge condition which, 
together with periodic boundary conditions, implies OqAo = 0. The same should be 
true, on average, when A is used. The behaviour of these two correlators is shown in 
Fig. | for the run W60b where the errors have been evaluated through jacknife. Fig. 
puts in evidence the flatness of (A A ) and (AqAo). As remarked above, (AoAo) has to 
be constant configuration by configuration, which is verified with a precision of ~ 10 -6 . 
More surprising is the fact that also (A A Q ) turns out to be constant configuration by 
configuration at the level of ~ 5%, also because in this case the behavior of the con- 
trol variable 9 , displayed in Fig. |l|, shows that A is far from satisfying the Lorentz 
condition on individual configurations. 

For both correlators, the error is just due to fluctuation of their constant (in t) 
value, configuration by configuration. 



In Table |l] we report the fit of the ratio: 

(■44) _ 

(AiAi) 



Cf(go) (17) 



as a constant in time. In Fig. || the Green function (.A-.A-) and the rescaled one 
Cf(go)(AiAi) are reported for the run W60b. The remarkable agreement between 
these two quantities confirms the proportionality shown in eq. (|i~0|). 

As shown in Table 0, Co(<?q) and Ci(go) coincide, within the errors, for the symmet- 
ric lattices W58 and W60a , while they have a different value for W60b and W64. This 
is probably due to the breaking of cubic symmetry. This interpretation is confirmed 
by considering the W60a (8 3 • 8) and W60b (8 3 • 16) lattices, with the same j3 value 
and configured so that the time extension of W60a is equal to the spatial extension of 
W60b. In fact we find that the coefficients Cb(<?o) evaluated from W60a and Ci(go) 
estimated from W60b agree within the errors. 
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The breaking of cubic symmetry, which depends on the ratio T/fV) 1 / 3 , could be a 
potential source of systematic error in the non-perturbative evaluation of renormaliza- 
tion constants on asymmetric lattices. 

The remaining Green functions, eqs.(p^-p^), which contain the insertion of d^A^ 
exhibit the expected behavior: the ones formed with A„ fluctuate around zero at a level 
of 10~ 5 , while those built through A'^ fluctuate around zero at a level of 1CT 2 - 10" 3 . 

We are now ready to show why the discrepancy between the values of 9, relevant 
to control the gauge-fixing algorithm, and the expectation values of 9 , is natural. In 
fact we have: 

9 = a 4 J d A x (d^(x)f (18) 
9' = a 4 J d 4 x (d^(x)f 

In other words, 6 and 9' are proportional to two different discretizations of the same 
continuum operator (d^A^) 2 . However, while 9 vanishes configuration by configuration, 
as a consequence of the gauge fixing procedure, 9' is proportional to (d^A'^) 2 , which 
has the vacuum quantum numbers and mixes with the identity. The expectation value 
of (d^A'^) 2 , therefore, diverges as so that 9 will stay finite, as a —* 0. 

4 Conclusions 

In this paper we investigated the problem of the definition of lattice operators con- 
verging, as a — » 0, to the fundamental continuum gauge fields. This construction is 
affected, at the regularized lattice level, by an enormous redundancy. However we 
found convincing non-perturbative evidence, based on numerical simulations, that this 
redundancy will be completely compensated, in the continuum, by the wave function 
renormalization needed in order to define finite Green's functions with gauge field in- 
sertions. Although, on general field theoretical grounds, the validity of such results is 
not unexpected, we stress the remarkable fact that it holds true also in this particular 
situation in which gauge-fixing is naively performed, disregarding the problems related 
to the existence of lattice and continuum Gribov copies. 

We have also discussed and solved the problem of the large discrepancy between the 
values of the control functionals 9 and 9 in terms of renormalization of power-divergent 
composite operators. 

A direction dependent effect in the value of the renormalization constants has been 
found on non-cubic lattices, which should caution against too naive application of 
infinite volume results to finite volume lattice numerical data. 

We believe that these features have a general validity and will survive a more 
thorough treatment of the gauge-fixing problem. 
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W58 W60a W60b W64 


# Confs 
Volume 


5.8 6.0 6.0 6.4 
20 100 50 30 

6 3 x 6 8 3 x 8 8 3 x 16 8 3 x 16 


Ci(g Q ) 
CoQfo) 


0.689(3) 0.729(1) 0.729(2) 0.757(2) 
0.690(7) 0.729(1) 0.750(1) 0.784(2) 


a" 1 


1.333(6) 1.94(5) 1.94(5) 3.62(4) 



Table 1: Summary of the lattice parameters used and relative values of Co and C, obtained as 
discussed in section ||. The a~ l values are taken from [19, |2(| and are computed through the string 
tension. 
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Figure 1: Typical behaviour of 6 and 9' vs gauge fixing sweeps at (3 = 6.0 for a thermalized SU(3) 
configuration 8 3 • 16. 
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Figure 2: Comparison of the matrix elements of (AoAo)(t) (open circles) and {A' A' Q )(t) (crosses) 
as function of time for a set of 50 thermalized SU (3) configurations at [3 = 6.0 with a volume 
V ■ T = 8 3 • 16 (run 1/F606); the errors are jacknife. 
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Figure 3: Comparison of the matrix elements of (A[A^)(i) (crosses) and the rescaled (AA) ■ Cf (<7o) 
(open circles) as function of time for a set of 50 thermalized SU (3) configurations at (3 = 6.0 with 
a volume V ■ T = 8 3 • 16 (run W60b). The data have been slightly displaced in t to help eye, the 
errors are jacknife. 
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